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ABSTRACT 


The magnetic activity of a star — which modulates the stellar wind outflow — shapes the immediate 
environments of orbiting planets and induces atmospheric loss thereby impacting their habitability. 
We perform a detailed parameter space study using three dimensional magnetohydrodynamic simu- 
lations to understand the effect of changing stellar wind magnetic field and planetary magnetic field 
strengths on planetary magnetospheric topology and atmospheric losses. It is observed that the rela- 
tive strengths of stellar and planetary magnetic fields play a significant role in determining the steady 
state magnetospheric configuration and atmospheric erosion. When the stellar field is strengthened or 
the planetary field is weakened, stellar magnetic field accumulation occurs at the day-side of the planet 
which forces the magnetopause to shift closer to its surface. The magnetotail opens up leading to the 
formation of Alfvén wings in the night-side wake region. We demonstrate how reconnection processes 
and wind conditions lead to the bifurcation of the magnetotail current sheet. With increasing stellar 
wind magnetic field strength, the day-side reconnection point approaches the planet thereby enhancing 
mass loss. We establish an analytic equation which successfully captures the modeled mass-loss rate 
variations of planets with changing magnetic field strengths. Our results are relevant for understand- 
ing how the interplay of stellar and planetary magnetism influence (exo)planetary environments and 
their habitability in star-planet systems with differing relative magnetic field strengths, or in a single 
star-planet system over the course of their evolution with age. 


Keywords: Unified Astronomy Tesaurus concepts: Stellar winds (1636); Magnetohydrodynamical sim- 
ulations (1966); Planetary magnetospheres (997); Exoplanet atmosphere (487); Exoplanet 
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1. INTRODUCTION 


A star interacts with the planets it hosts via magnetized stellar wind which directly impacts their magnetospheric 
configuration and atmospheric loss. Gravitational (tidal) interactions can be neglected if the planet lies far away (as in 
the case for Sun-Earth system) and in such a scenario, the stellar wind independently affects the planetary environment 
without the possibility of the coronal structure of the star getting modified. The evolution of planetary atmospheres 
and the nature of its interaction with the stellar wind are profoundly affected by changes in their respective magnetic 
field strengths (Cohen et al. 2015; Basak & Nandy 2021). Discovery of several exoplanets (Mayor & Queloz 1995; 
Pepe et al. 2014; Lunine et al. 2009) has sparked interest in exploring the signatures of life-sustaining conditions and 
understanding habitability from the perspective of the origin and evolution of planetary atmospheres (Lammer et al. 
2009; Lammer 2013; Pollack & Yung 1980; Cridland et al. 2017). Habitability — in the astrophysical context — depends 
on the ability of the stellar wind forced planet to hold on to its atmosphere (Alvarado-Gémez et al. 2020) and some 
other aspects; this work focuses on stellar wind planetary interactions. 

The age of the star determines its activity which dictates the properties of the outflowing wind. Observations confirm 
the varying magnetic activity of the Sun during its evolutionary phases (Nandy & Martens 2007; Vidotto 2021), as 
well as variations in field strength of solar system planets and stars outside our solar system (Stevenson 1983; Kiefer 
et al. 2017; J.-D. do Nascimento et al. 2016). Exploration of magnetic fields on exoplanets is a research pursuit of 
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high topical interest (Oklopéié et al. 2020; GrieBmeier 2015). With the detection of diverse magnetic activity in a 
number of star-planet systems, the question that naturally arises is — what is the consequence of magnetic interactions 
of (exo)planets with their host star? This is a critical question because the answer to this determines whether a planet 
can retain its atmosphere during its long-term interaction with the harboring star and consequently, its habitability 
time-scale. We explore how the nature of this interaction is determined by the relative strengths of the stellar wind 
magnetic field and planetary magnetosphere. Several efforts have been made to understand how stellar winds influence 
the planetary magnetosphere as well as its atmosphere (Lammer et al. 2012; See et al. 2014; Vidotto & Cleary 2020; 
Harbach et al. 2021; Nandy et al. 2021). Variation in stellar magnetic activity (Nandy 2004; Nandy & Martens 2007; 
Tripathi et al. 2021; Brun et al. 2014; Vidotto 2021; Nandy et al. 2021) introduces variations in stellar radiation (Spina 
et al. 2020), stellar wind speed (Finley et al. 2018), magnetic field strength of plasma winds (Vidotto et al. 2015) and 
may lead to magnetic storms as well. As a consequence of these interconnected phenomena, the planetary dipolar field 
is deformed resulting in a magnetospheric structure that may vary from what we often observe on the Earth (Gallet 
et al. 2016). 

In this paper, we perform a thorough parameter space study using numerical modelling for understanding interactions 
in star-planet systems and the effect of stellar activity evolution on planets with different magnetospheric strengths. 
The magnetized stellar wind and the planet’s inherent magnetic field are projected to have a significant impact on 
planetary habitability and atmospheric mass loss (Nandy 2004; Khodachenko et al. 2008; Nandy et al. 2017) . We carry 
out three dimensional global magnetohydrodynamic (MHD) simulations using a broad range of stellar and planetary 
magnetic field values for providing a comprehensive picture of the interaction process. Most extrasolar giant planets 
that have been discovered so far are anticipated to have significant ionospheres. To simplify the global simulation 
and based on prior studies (Koskinen et al. 2010; Yan et al. 2019; Gronoff et al. 2011), we consider the planet to be 
surrounded by a perfectly conducting plasma atmosphere. We study the impact of increasing stellar wind magnetic 
field on magnetopause stand-off distance and explore the conditions that lead to the magnetic pile-up in the day-side 
region. Another significant occurrence we examine is the formation of Alfvén wings where the magnetotail opens up 
when the upstream wind Alfvénic Mach number is low and the current sheet length shortens and bifurcates in the 
night-side region of the planet. We find a criterion for current sheet bifurcation that is consistent with all cases of 
our study and establish an analytical relation for the variation of mass loss rate with changing stellar and planetary 
magnetic fields. 

This study is relevant for any star-planet system and in those systems in which the magnetic fields of either or 
both entities vary considerably with time. Our findings are significant for a better understanding of (exo)planetary 
atmospheres and determining the impact of varying magnetic field strength on habitability. 

This paper begins with the model description in section 2 in which we give an overview of the theory and numerical 
setup employed in this study. The detailed findings are presented in section 3 followed by the conclusion in section 4. 


2. MODEL DESCRIPTION 


We adapt the Star-Planet Interaction Module (CESSI-SPIM) developed by Das et al. (2019) for simulating different 
configurations of stellar wind and planetary magnetospheres. The governing set of resistive MHD equations are given 


by: 
O.p +V - (pv) =0 (1) 
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where the symbols p, v, B, J, P, E, and g denote density, velocity, magnetic field, current density, pressure, total 
energy density and gravitational acceleration due to the planet respectively. The expression for total energy density 
is given by 


E= | | (5) 


for an ideal gas equation of state. 

The computational domain extends from -80 Rp to 200 Rp in the x-direction, -45 Rp to 45 Rp in the y-direction 
and -200 R, to 200 R, in the z-direction, where R, is the radius of the planet which is located at the origin of the 
Cartesian box. The region extending from -2 Rp to 2 Rp is resolved by 12 grids in all three directions i.e. 1 Rp is 
resolved using 3 grids in the planetary vicinity. Keeping in mind our modest computational facility, a combination of 
stretched and uniform grid types is used in the x and z directions. In the x-direction, the regions extending from -80 
Rp to -2 Rp and 2 Rp to 10 R, are resolved using 156 and 16 grids respectively i.e. 1 Rp is resolved using 2 grids. 
Grids with stretching ratio 1.005175 are used from 10 Rp to 200 Rp. In the y-direction, the regions from -45 Rp to 
-2 Rp and from 2 R, to 45 Rp are resolved by 86 grids each i.e. 1 Rp is resolved using 2 grids. In the z-direction, 
the regions from -22 R, to -2 Rp and 2 R, to 22 R, are resolved by 40 grids each i.e. 1 Rp is resolved using 2 grids 
whereas grids with stretching ratio 1.004958 are used for the regions extending from -200 R, to -22 R, and 22 R, to 
200 Ry. 

In this study, an Earth-like planet is considered with similar mass, radius and tilt while the stellar wind is assumed 
to have a speed of 270 km/s corresponding to quiet times. We perform a set of simulations by varying the magnetic 
field strength of both the planetary magnetosphere (0.1Be,0.5Be, 1Be, 2Be and 5Be) and stellar wind (0.1 nT, 0.5 nT, 
1 nT, 2nT, 5 nT, 10 nT, 30 nT, 50 nT and 75 nT) with all possible combinations. Here Be = 3.1 x 104 nT denotes 
surface equatorial magnetic field strength of the Earth’s dipolar field. 


Table 1. Values of physical parameters used in the simulations and their respective notations. 


Physical quantity Notation Value used 
Density in ambient medium Pamb 1.5 x107 gcm’ 
Pressure in ambient medium Pamo 2.49 x10771} dyne cm~? 
Density of stellar wind Psw 4 Pamb 
Velocity of stellar wind Vsw 2.7 x10’ cms + 
Adiabatic index y 5/3 
Planetary mass My 5.972 x10?” g 
Planetary radius Rp 6.371 x10 cm 
Intrinsic planetary magnetic field Bpi 0.1Be - 5Be * 
Stellar wind magnetic field Bai 0.1 aT - 75 nT 
Magnetic diffusivity n 1073 cm? s7! 
Magnetospheric tilt angle Opi 11° 


* The symbol Be represents the magnetic field for the case of the Earth’s dipole. 


In order to study the effect of stellar activity evolution on planetary atmospheric loss which directly impacts hab- 
itability, we initialize a conducting plasma atmosphere surrounding the planet. The density profile in the vicinity of 
the planet is defined by 

Ppl = 1O Pair r < Rp 


parm(r) = pp + CEP [rannd o( 2 —2)} +1) 


Rp <r < Sites 


(6) 


where pp; and Pam» are densities of the planet and ambient medium respectively and r is the radial distance from the 
origin. For more details on the justifications for the choice of the above atmospheric profile, interested readers may 
refer to Sec. 2.1 of Basak & Nandy (2021). The pressure distribution in the atmosphere is evaluated by numerical 


integration of the equation 
dP 
— = ~parm(r)o(r) (7) 


Here, the gravitational field is given by g(r) = — CNet where Mp is the planetary mass. The pressure inside the planet 


is evaluated by extrapolating the value of pressure at planet-ionosphere boundary. The density and pressure in the 
region r > 3Ry is initialized to be equal to that in the ambient medium (Table 1). 


The input parameters at the stellar wind injection boundary (at x= -80 Rp in yz plane) are obtained by solving the 
Rankine-Hugoniot magnetized jump conditions for the given shock velocity. A southward oriented (SIMF) stellar wind 
is injected with density Psw = 4 Pamp. For all other boundary faces of the Cartesian box, force-free outflow boundary 
condition is implemented. The magnetic fields of the wind and (or) planet are varied in each run and the steady state 
configuration of the planetary magnetosphere and atmospheric loss are analyzed. 


3. RESULTS AND DISCUSSION 


When the stellar wind impinges upon the planetary magnetosphere, a dynamic pressure balance is achieved at the 
day-side of the planet resulting in the formation of a magnetopause which inhibits the penetration of stellar plasma 
and protects the inner atmosphere from erosion. The balance between thermal, magnetic and dynamic pressures 
determines the shape and location of the magnetopause. An earlier study by Basak & Nandy (2021) considers the 
magnetopause stand-off distance (Rm ,) to be the location of balance between the dynamic P4 (incoming wind) and 
magnetic pressure P, (near the planet) curves. This approach works well if the stellar wind magnetic field (Bsw) 
is weak since the magnetic pressure of the incoming wind is very small and may be neglected in the calculations. 
However, when the stellar magnetic field is relatively stronger, it accumulates outside the magnetopause and the 
magnetic pressure in the magnetosheath region cannot be ignored any more. Our study shows that the point of 
intersection between Pm and Py migrates outward nearer to the bowshock as B,y is increased (Figure 1). For a 
strongly magnetized wind, it is found that the dynamic pressure Py is sufficiently high even after this intersection 
point which indicates that stellar wind plasma is able to penetrate beyond this location and therefore, it cannot be 
regarded as the magnetopause stand-off distance. Thus, we need a different approach for evaluating the magnetopause 
stand-off distance for the case of strongly magnetized winds. It has been shown in previous studies (Schield 1969; 
Shue & Chao 2013; Lu et al. 2015) that at the magnetopause, the sum of magnetic and thermal pressures (denoted 
by Pint+m henceforth) on the planetary side balances the incoming stellar wind (see also Shue & Chao 2013). These 
considerations result in the relation 


[Pan + Pm] = [Pant Pr! (8) 


planet—side star—side . 
The distance from the planet where the Pih+m curve peaks and its gradient becomes zero along the subsolar line is 
then considered to be the magnetopause stand-off distance. 

The intersection of P;, and P4 curves gives the location of bow shock stand-off distance. Figure 1 shows the pressure 
balance for obtaining bow shock distance (Rs) and the magnetopause stand-off distance (Rmp) for the case of planetary 
magnetosphere Bp = Be (Be = 3.1 x 10* nT) with varying stellar wind magnetic field. Note that only southward 
interplanetary magnetic field (SIMF) is considered for this study. The total thermal pressure (P;;,) are estimated from 


— Pswt 


the input parameters (Mignone et al. 2007) while dynamic pressure of stellar wind (P4 = Psuv?) and magnetic pressure 


Tu= 2) are evaluated from evolving simulation variables. As the stellar wind magnetic field strength is increased, 
the magnetopause moves closer to the planet and the bow shock moves further away. For weakly magnetized stellar 
wind as shown in Figure 1(a)-(c), the stellar field accumulation is not significant and the intersection of P4 and Pm 
gives the Rmp location as depicted in Basak & Nandy (2021). As Bs y increases further, it is evident that stellar 
wind magnetic pressure builds up significantly in front of the planet (Figure 1 (d)-(f)). With the increase in Bsw, 
both thermal and magnetic pressure increase in the magnetosheath region initially up to some distance. However, the 
enhancement in magnetic pressure gets prominent gradually due to the accumulation of stellar wind magnetic field 
outside the magnetopause. 

Just after the magnetopause, magnetic pressure drops to nearly zero representing the day-side reconnection X point. 
Thermal pressure rapidly increases near the magnetopause due to energy conversion in the magnetic reconnection 
region. (see e.g. Lu et al. 2015) 

Pressure balance for different planetary magnetic field strength (B,) keeping the stellar wind magnetic field constant 
reveals that with decreasing B, both bow shock and magnetopause move towards the planet (Basak & Nandy 2021). 
This implies that for lower values of B,, stellar wind is able to penetrate closer to the planet resulting in greater 
atmospheric loss. We find significant increment of magnetic pressure in the magnetosheath region for the case of 
weaker planetary magnetic field as compared to stronger field cases. Thus, either strengthening the stellar wind 
magnetic field or weakening the planetary magnetic field leads to magnetic field accumulation at the day side of the 
planet. Hence, the ratio of the intrinsic planetary magnetic field and stellar wind magnetic fields is an essential factor 
which determines the dynamics of the star-planet interaction. This is in keeping with expectations. 
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The interaction of plasma winds and planetary magnetosphere results in certain deformation in the planetary dipolar 
field. This deformation leads to a tear-drop shaped magnetospheric structure that we usually observe for the Earth. 
A series of simulation results of interaction of planetary magnetosphere with stellar wind magnetic field is shown in 
Figure 2. In each subplot, magnetic field lines (white sreamlines) are traced over the current density (colormap) in the 
y = 0 plane with stellar wind plasma flowing from left to right. In plots (a)-(c), the simulations show the formation of 
extended magnetotail structures for varying planetary magnetic field Bp = 0.1 Be, 1 Be, and 5 Be respectively. Stellar 
wind magnetic field strength is fixed at Bsw = 2 nT for all cases. As Bsw is increased to 10 nT [panels (d)-(f)] for 
the same set of planetary magnetic field values, the magnetotail begins to open up and on further increase of Bsw to 
50 nT [panels (g)-(i)] the lobes of the tail open up even further (Turnpenney et al. 2020). The Alfvén Mach number 
is an important control parameter that governs the behaviour of plasma embedded in magnetic field. With increasing 
stellar wind magnetic field strength, the upstream Alfvén Mach number decreases which leads to the appearance 
of Alfvén wings across the planet (Ridley 2007). Figure 3 illustrates the sub-Alfvénic plasma interaction with the 
magnetized planet and the formation of three dimensional Alfvén wings in the planetary magnetotail as shown by the 
region enclosed by Alfvénic Mach number M4 isosurface (M4 < 0.3). When stellar wind magnetic field is strong, 
field lines bend but do not strongly drape the planet. Chané et al. (2012) have shown the observational evidence 
of Alfvén wings formation on the Earth during 24 and 25 May 2002. Several moons in our solar system including 
Ganymede, Io, Europa, and Enceladus have been observed to form Alfvén wings. Saur, J. et al. (2013) have found 
signatures of sub-Alfvénic plasma interactions in several star-(exo)planet systems. With decreasing upstream Alfvén 
Mach number, wings form at larger angles from the equatorial plane and the length of the current sheet becomes 
shorter i.e., reconnection takes place closer to the planet in the magnetotail region as shown in column-wise panels of 
Figure 2. The row-wise panels of Figure 2 show that for a given stellar wind magnetic field strength, Alfvén wings are 
more prevalent in case of weaker planetary magnetic field. Here, we present only a few simulation results to show how 
the magnetotail opens up. A similar pattern is found for other simulations that involve all possible combinations of 
parameters (described in Table 1). Vernisse et al. (2017) have shown that interactions which involve either a weakly 
magnetized obstacle or sub-Alfvénic upstream plasma velocities or both can lead to the dominance of an Alfvén wing 
structure — and this is corroborated by our simulations. 

To analyze the opening up of the magnetotail in greater detail, we locate the distance from the planet where the 
current sheet bifurcates in the night-side region. Figure 4(a) shows the current density plot in xz-plane (left) and 
corresponding current density magnitude (Jmag) variation along the subsolar line (right). Since the current sheet has 
a finite thickness along y and z directions, we consider a hypothetical rectangular box outside the planet’s atmosphere 
encapsulating the current sheet and extending from -0.5 R, to 0.5 R, along both y and z directions. Thereafter, the 
spatially averaged Jmag is calculated on yz slices of the cube and its variation along x-direction is plotted. We find 
that the current sheet starts to bifurcate typically when the current density magnitude drops off to about 45% - 50% of 
its peak value. These are represented by line B and line A respectively in Figure 4(a). This result is robust considering 
all different sets (Table 1) of our simulations. The current sheet, extending from the planet to its bifurcation point on 
the night-side region, is hereby referred to as the current sheet length. Variation of current sheet length with planetary 
& stellar wind magnetic field is shown in Figure 4(b). As the planetary magnetic field B, increases, the bifurcation 
point recedes away from the planet. For Bp = 2B. we do not see any bifurcation within the simulation domain until 
Bsw = 30 nT. For a stronger planetary magnetic field of Bp = 5Be, bifurcation starts around Bsw= 50 nT onwards. 

It is evident from Figure 4(b) that the current sheet bifurcation point shifts nearer to the planet when either the 
stellar wind magnetic field increases or the planetary magnetic field weakens. 

We also investigate how the variations in stellar wind and planetary magnetic field strengths impact the mass- 
loss rate of planetary atmospheres. To compute the total mass-loss rate, we consider a cube of edge length 6.6 R, 
(extending from -3.3 R, to 3.3 Rp in all three directions) with its origin coinciding with the centre of the planet such 
that the whole planet including its atmosphere is enclosed within the cube. Mass loss rates are computed across all 
six faces of the cube which contribute to the total mass loss. The rate of atmospheric mass loss depends on the extent 
of stellar wind penetration into the planetary atmosphere. This is essentially denoted by a point at the dayside — 
the reconnection X point — beyond which planetary magnetic field lines merge with the stellar field. If the magnetic 
field strength of the stellar wind is increased gradually keeping that of the planet fixed, the X-point is expected to 
shift closer to the planet inducing greater atmospheric loss. However, the distance of the X-point from the planet 
does not decrease continually with increasing stellar field. The reconnection neutral point saturates at a certain level 
and then onwards, the stellar field starts accumulating at the day-side [Figure 1 (d)-(f)] similar to that of an imposed 
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magnetosphere; at this stage mass loss rate attains a limiting value. Figure 5 (a) confirms increase in atmospheric 
mass-loss rate with increasing stellar wind magnetic field strength. Figure 5 (b) shows the plot of the X-point distance 
from the planetary surface with increasing stellar wind magnetic field for different values of the planetary field. This 
signifies that for a weak planetary field, the X-point distance saturates to a steady value at lower stellar magnetic 
fields while for a stronger planetary field, the saturation occurs for higher values of the stellar wind field. Figure 5 
(a) shows that for a weaker planetary field, the loss rate increases initially and then saturates at lower values of the 
stellar wind field while for the stronger planetary field, the loss rate requires higher values of stellar wind field to attain 
saturation. The atmospheric loss rate is found to have a strong negative correlation with the X-point distance from 
the planet. The Pearson correlation coefficient values are -0.9800, -0.9683, -0.9753 and -0.9659 for planetary magnetic 
fields 0.5B., 1B., 2Be, and 5B, respectively. 

Figure 5(a) also indicates that for a conducting plasma atmosphere, as considered in the present study, the mass 
loss rate increases with increasing strength of the planetary magnetosphere. This is because a stronger magnetosphere 
provides a larger interaction area with the ambient medium and therefore, may lead to greater loss when it interacts 
with the stellar wind. Moreover, as the atmospheric matter in our study is not considered to be charge neutral, it 
is expected that the atmospheric plasma will take part in magnetic reconnections and be lost in the process. Earlier 
studies show that the atmospheric matter stretches out and escapes through polar caps and cusps during reconnections 
(Gunell et al. 2018; Egan et al. 2019; Sakata et al. 2020; Carolan et al. 2021) resulting in a larger atmospheric mass 
loss for a stronger magnetosphere. Gunell et al. (2018) shows that an intrinsic magnetic field does not necessarily 
protect a planet against atmospheric escape and the escape rate can be higher even for a highly magnetized planet. 
Therefore the role of intrinsic magnetic field in the protection of planetary atmosphere is a complex question. 


SEMI-ANALYTICAL EXPRESSION FOR PLANETARY MASS LOSS RATES 


Our study shows that the global dynamics of the planet, current sheet length and atmospheric mass loss rates are 
significantly influenced by the relative strength of the planetary and stellar wind magnetic fields. For different sets 
(Table 1) of our simulations, Figure 6 shows that ratio of planetary to stellar wind magnetic fields is an important factor 
in governing atmospheric mass loss rates. Based on theoretical considerations we establish an analytical relationship 
between these two in order to ascertain the dependence of mass loss rate on the variation of planetary and stellar wind 
magnetic fields. 

The unpeturbed planetary dipolar magnetic field falls as 4, thus pressure balance at magnetopause can be written 
as 


B? k?, B2 RS 
Ddynsw t+ + Pthsw = ~~ + inp; (9) 
210 2uoR$,p 


where dynamic pressure due to stellar wind, thermal pressure due to stellar wind and thermal pressure due to planetary 
magnetosphere at magnetopause are depicted by Pdynsw and Pthsw, Pthp respectively. The symbols Bsw and Bp 
represent the stellar wind and planetary magnetic field strengths respectively while Rmp is the magnetopause standoff 
distance in terms of planetary radius Rp. The factor km is a measure of the compression of the dipolar magnetic field 
by the stellar wind at the magnetopause (Mead & Beard 1964; Tsyganenko 2005). 

From equation 9, the magnetopause stand-off distance can be expressed as 


(10) 


where C = 2uo(Pdyn sw + Pthsw — Pthp). As pointed out in section 3, the atmospheric mass loss rate (M ) is linearly 
anti-correlated with the day side X-point distance (Xpoint) and can be expressed as M = —AX point + B where A 
and B are fitting parameters that need to be determined. Moreover, in case of an SIMF the magnetospheric stand-off 
distance, Rmp is approximately equal to the day side reconnection point (Xpoint). Therefore, equation 10 can be 
expressed as 


-B (11) 


T 


Equation 11 gives an analytical relationship between mass loss rate and the ratio of planetary and stellar wind magnetic 
field strengths. 

Figure 6 shows that our modelled output of atmospheric mass loss rate fits reasonably well with the analytical 
expression given by 11. The value of Rsquare (goodness of fit) for the modelled fit is greater than 0.99 for all instances 
of planetary magnetic fields. 


4. CONCLUSIONS 


In this study, we have used 3D global magnetohydrodynamic modelling to explore the effect of varying stellar 
and planetary magnetic field strengths on the steady state magnetospheric configuration and atmospheric loss rates of 
planets. This study is important for understanding how the magnetic activity evolution of stars affects the environment 
and habitability of planets with different magnetospheric strengths. Our results are therefore applicable to a wide 
domain of far-out (exo)planetary systems. 

Computing pressure balance at the day-side of the planet illustrates that either strengthening the stellar wind 
magnetic field or weakening the planetary magnetosphere results in stellar field accumulation in front of the planet 
similar to that of an imposed magnetosphere. In such a scenario, the magnetopause stand-off distance cannot be 
obtained using the conventional procedure of balancing the dynamic pressure of the incoming wind and magnetic 
pressure close to the planet but rather the sum of both magnetic and thermal pressures need to be considered. It is 
found that for a particular strength of planetary magnetopshere, increasing the magnetic field of stellar wind makes 
the bowshock shift away from the planet while the magnetopause is formed closer, implying greater wind penetration. 

We find that the relative strength of planetary and stellar wind magnetic fields plays a critical role in determining 
the magnetic reconnection in the magnetotail region. The long extended planetary magnetotail which exists for a 
moderately magnetized stellar wind or strong intrinsic planetary magnetic field ceases to exist for cases when the 
stellar wind magnetic field becomes extremely strong (i.e. the Alfvénic Mach number decreases) or the planetary 
magnetic field becomes very weak. This leads to the formation of Alfvén wings in the night-side wake region. Our 
findings suggest that the magnetotail begins to open up and the point of bifurcation shifts closer to the planet when 
either the magnetic strength of the stellar wind is increased or the intrinsic planetary magnetic field is reduced. We 
observe that the magnetotail current sheet starts to bifurcate when current density magnitude drops to about 45-50 
% of its maximum value. 

We also find that the atmospheric mass-loss rate increases with increase in stellar magnetic field due to greater wind 
penetration close to the planet. The atmospheric loss rate saturates at lower values of stellar wind magnetic field for 
a weak planetary field. In the case of a stronger planetary field, however, a stronger stellar wind magnetic field is 
required for the saturation of the atmospheric loss rate. The variation of the day-side X-point distance from the planet 
with stellar wind magnetic field provides a thorough explanation for the saturation of atmospheric loss rate since the 
day side X-point distance from the planet and atmospheric loss rate are found to have a strong negative correlation. 
By modelling the impact of planetary magnetic field on atmospheric escape processes, we corroborate the study by 
Gunell et al. (2018) that the escape rate can be higher for strongly magnetized planets. 

We establish an analytical relationship between the mass loss rate and the ratio of planetary and stellar wind 
magnetic fields. Our study shows that the global dynamics of the planet, current sheet length, and atmospheric mass 
loss rate are significantly influenced by the relative magnetic field strengths of the star and the planet. The results of 
our numerical simulation indicate that the analytical expression derived in this study describes quite well the modeled 
atmospheric mass loss rates. 

The detailed parameter space study presented in this paper will be helpful for identifying stellar and planetary 
conditions that lead to a particular magnetospheric configuration in solar and exoplanetary systems. The results also 
illustrate the impact of magnetic field variability on atmospheric loss rates which is relevant for understanding the 
habitability of (exo)planets. 
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A. APPENDIX INFORMATION 


(b) BSw = 2nT (c) Bsw = 5nT 


fe) 

A 

Cc 

sat 

sa 

(0.0) 

x -40 -30 -20 -10 0 -40 -30 -20 -10 0 -40 -30 -20 -10 0 
~ (d) Bsw = 10n (e) Bsw = 30nT (f) Bsw = 50nT 
@M 10 10 10 

km 

mD) 

Ww 

WY 

cD) 

_ 

ial 


Figure 1. Pressure balance for bowshock and magnetopause: Plots of different pressures in the planetary day-side for magne- 
tospheric field B, = Be (3.1 x 104 nT) with varying stellar wind field strength (Bs). The symbols Pin, Pm, Pa, Pin+m denote 
thermal, magnetic, dynamic and sum of thermal and magnetic pressures respectively. The symbols Ry, and Rmp denote the 
positions of the bowshock and the magnetopause stand-off respectively. The value of R,, is obtained by balancing P4 (wind) 
and Pin (magnetosheath) while Rmp is evaluated by identifying the location close to the planet where the gradient of Pinim is 
Zero. 
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(a) Bp = 0.1 Be, Bsw = 2nT (b) Bp = Be, Bsw = 2nT (c) Bp = 5 Be, Bsw = 2nT 
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Figure 2. Steady state planetary magnetospheric configurations for varying stellar and planetary magnetic fields. The back- 
ground colormap shows current density while magnetic field is plotted using white streamlines in the y = 0 plane. The planet 
is located at the origin and the stellar wind flows from left to right. Stellar wind magnetic field increases along each row (top 
to bottom) while planetary magnetic field strength increases along each column (left to right). Distances are measured in units 


of the planetary radius. 
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Figure 3. Steady-state magnetospheric configuration of an (exo)planet with Bp = 0.1B. and Bsw = 50 nT. The white surface 
encloses regions where Mach number Ma < 0.3. Pink and yellow streamlines represent magnetic field lines originating from the 
planet and pure stellar wind respectively. 
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Figure 4. Current sheet bifurcation due to magnetotail opening: (a) Plot of current density in the y = 0 plane (left) and 
corresponding variation of its magnitude (Jmag) along the subsolar line (right) for Bp = Be and Bsw = 30 nT. Line A denotes 
the location of maximum Jmag while line B marks the bifurcation point of the current sheet. (b) Left: Variation of current 
sheet length with stellar wind magnetic field for three different cases of planetary magnetic field strength. Right: Illustration of 
current sheet length variation with stellar and planetary magnetic field strength. The light grey region represents cases where 
the current sheet length has not been bifurcated within the simulation domain. 
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Figure 5. Dynamics of planetary atmospheric loss: Plots of (a) total mass loss rate from a cube of side 6.6 Rp with the planet 


at its centre and (b) 


variation of day-side X-point distance from the planetary surface with stellar wind magnetic field strength 


for different magnetospheric strengths. 
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Figure 6. Dependence of atmospheric mass loss rate on the ratio of the planetary to the stellar wind magnetic field strength. 


Simulation results are represented by square boxes, while the fit from the analytic expression (equation 11) is shown by dotted 
lines. 


